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ABSTRACT 






CNJ I Aims. For the first time, the lunar occupation technique has been employed on a very large telescope in the near-IR with the aim of 

achieving systematically milliarcsecond resolution on stellar sources. 

Methods. We have demonstrated the burst mode of the ISAAC instrument, using a fast read-out on a small area of the detector to 
record many tens of seconds of data at a time on fields of few squared arcseconds. We have used the opportunity to record a large 
number of LO events during a passage of the Moon close to the Galactic Center in March 2006. We have developed and employed 
Q , for the first time a data pipeline for the treatment of LO data, including the automated estimation of the main data analysis parameters 

5^ . using a wavelet-based method, and the preliminary fitting and plotting of all light curves. 

"^ ' Results. We recorded 5 1 LO events over about four hours. Of these, 30 resulted of sufficient quality to enable a detailed fitting. We 

jrt ' detected two binaries with subarcsecond projected separation and three stars with a marginally resolved angular diameter of about 

2 milliarcseconds. Two more stars, which are cross-identified with SiO maser, were found to be resolved and in one case we could 

recover the brightness profile of the extended emission, which is well consistent with an optically thin shell. The remaining unresolved 

Cn ' stars were used to characterize the performance of the method. 

^ ' Conclusions. The LO technique at a very large telescope is a powerful and efficient method to achieve angular resolution, sensitivity, 

^O , and dynamic range that are among the best possible today with any technique. The selection of targets is naturally limited and LOs 

^^ . are fixed-time events, however each observation requires only a few minutes including overheads. As such, LOs are ideally suited to 

^^D ' fill small gaps of idle time between standard observations. 

CN ■ 

• ' Key words. Techniques: high angular resolution ~ Astrometry - Occultations - Binaries - Masers 

o 

oo , 

(*— ^ , 1. Introduction for observation (dominated by overheads, since an occultation 

• • . _, , , ,. , . , , . ,T ^N ,-, , , l^sts much less than 1 s) and for data analysis is significantly 

> The method of observing unar occultations (L.O) of background ^^^^^^^ ^^^^ ^^^ ^^^^ ^^^^ ^.^ ^^^ j^^ resolution methods. 

r'^ stars to derive their angular diameter, as well as other charac- _j, r^n-j-i .c-^z i 
S^ • ^ ■ ,. ui--fui- ijrjj J Or course, LO suiter from a number or significant draw- 
rN teristics such as binanty, has been employed tor decades and , , jz , , , , ■ 
Lj ' . , , , ,. , ,, •,,• J / X 1 1 backs, first among them that the sources cannot be chosen at 
13 : provides angular resolution at the milliarcsecond (mas) level •„ 4-, ., , , ir.^ .c . i ■ i i 
Cu r , J ,.? r .u J- . r .u . 1 J T-u- will. The Moon covers only about 10% ot the celestial sphere 
■ - - independently ot the diameter ot the telescope used. This sur- . . , • t ^ jz i • i i i 
.u ji .• T -^ c .u 1 ; • 1 . 1 in Its apparent orbit. LO are fixed time events, and as such sub- 
passes the dinraction limit ot even the largest single telescopes . , , . , , ■ V- n i i 

, . , ^, , . „, , T ■ ^ c . /T nT\ lect to weather and instrumental downtimes. Finally, the lunar 

and rivals the resolution ot long -baseline interterometry (LBI) ,., , ., ,„ „, iii 

.^, , ,. c T. J J c . n • 11 li, T /-> limb only provides a 1-D scan ot the source, although observa- 

even with baselines ot hundreds ot meters. Basically, the LO . „ -^ ^,.„ . ,.„ , .... 

^ , . ,. ^, , ,. , j.rr. ,. J ,, tions from diiterent sites, or at ditterent epochs, can in principle 

technique relies on the lunar limb as a diitracting edge, rather , ijji^ii- i^ i ■ 

^ ° ° be combined under favorable circumstances. Due mainly to the 



than on the diameter of the telescope. The diffraction fringes are 
generated in space, and due to their relatively fast motion over 
the telescope they must be sampled at rates of order 0.1-1 kHz. 

Combined together, these factors greatly diminish the adverse ,, ,. ,, , ,„ jr.ii- 

„, c . u-.i-i I.- t .u i-.r mance, the data analysis and the results can be found in Richichi 

efrects ot atmospheric turbulence, which are the main limit ot i,„„, | i ,„„,i , „■' i-^„,i _, ^,,.r..,^ , , 

,u u- u i w . u • r /-, 1 •. . J 1985 19961 ) and Fors (I2006I I. The CHARM2 catalogue alone 

other high angular resolution techniques. LO also permit not ^ „ . , ',. — , ! -„„-■ ,. — TTT, ^ t ^ i i ^ ,j r , ■ , 

. jij ..•:• 1 (Richichi et al. I2005I I hsts 1815 LO results in the field of high 



chromatic properties of the scattered light background around 
the Moon and of the diffraction fringes, the near-IR is ideally 
suited to observe LO. More details on the method, its perfor- 



angular resolution. 



only a model-dependent derivation of source parameters such 

as the angular diameter and binary parameters, but also a model- 

independent reconstruction of the brightness profile of the source ^hile the angular resolution of LO is set mainly by the lu- 

according to maximum-likelihood, or even a unique reconstruc- "^i" li"^^ '"^*^i" *^" ^y telescope size, this latter parameter ob- 

tion by light curve inversion in special cases. The time required ^lo^^^y ^as a crucial role in determining the limiting sensitivity. 

Richichi (1997 ) studied the performance of LO under a number 

Send offprint requests to: A. Richichi of circumstances, including the then-futuristic use of IR array 
* Based on observations made with ESO telescopes at Paranal detectors on a 8-m c lass tele scope. This ha s recently be come re- 
Observatory ality: Richichi et al. ( |2006bl l and Fors et al. (120061 |2008] l reported 
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Fig. 2. Color-magnitude diagram for the sources of which we 
recorded LO events as listed in Table[T](open squares). The small 
dots mark all sources occulted by the Moon in the same time 
span, to the limit K< 7.5 mag. 



in a preliminary form on the use of the ISAAC instrument in the 
so-called burst mode at the ESO VLT 8.2 m Antu telescope. 

In the present paper we provide a detailed account of the ob- 
servations carried out during a few hours in the night of March 
21, 2006, when 51 LO events were recorded. A second batch of 
observations was carried out in August 2006: due to their large 
number and their different nature, these sources will be discussed 
in a separate paper. In Sect.|2]we describe the observations and 
we provide details of the data reduction. In Sect.[3]we discuss the 
results, which include new binaries, resolved stars and the near- 
IR counterparts of two radio maser sources. We also character- 
ize the performance of this specialized observing mode, which 
combines the highest angular resolution and sensitivity possible 
today in the near-IR. In the conclusions we describe the integra- 
tion of the LO technique in the service mode operations scheme 
in place at the Paranal observatory and our plans for routine LO 
observations at that site. 



2. Observations and data analysis 

We observed a passage of the Moon in a crowded region of the 
Galactic Bulge, in the night of March 21-22, 2006. The cen- 
ter of this region was located at ^ 17H2™ and -28°29'. Fig.[T] 
shows the area, with the apparent lunar path superimposed. It 
can be seen that, at closest approach, the limb of the Moon as 
seen from Paranal reached only a few arcseconds from the true 
Galactic Center. The region is heavily reddened by interstellar 
dust. A near-IR color-magnitude diagram is shown in Fig. |2] 
Due to extinction, few sources in this region of the sky have op- 
tical counterparts. Our LO predictions were based mostly on the 
2MASS Catalogue (Cutii et al l2003T l. To the Hmit K< 7.5 mag a 
total of 509 sources were predicted to be occulted by the Moon 
in 5 hours. Our observations began at 6:11 UT and ended at 
10:25 UT. A total of 51 events were recorded. The efficiency was 
limited in the first part of the observations by several missed or 
wrongly recorded events. In fact, this run represented a commis- 
sioning of the burst mode at ISAAC, which at the time had not 
yet been extensively tested. In the last part of the observations, 
with the technical aspects under control, the typical interval be- 
tween recorded events was as short as three minutes and limited 
by the telescope pointing and data storing times. 



The limit K< 7.5 mag was quite arbitrary. In fact, in the same 
time span the number of occulted sources with K$ 1 1 mag was 
more than 15,000. As we will show later this magnitude limit 
would have been realistic, but it was clearly impossible to ob- 
serve hundreds, not to mention thousands, of sources. Even a 
manual selection of the best ones would have been a daunting 
task, and for this we developed a prioritization rule which takes 
into account the magnitude, the color (redder objects were given 
more weight), and the time intervals between the events. We also 
flagged and gave additional priority to sources with previous ob- 
servations or classification avalailable in the SIMBAD database. 

The lunar phase was -54%. In general LO observations are 
possible only on the dark lunar limb, thus we observed reap- 
pearances: the telescope was pointed at the nominal star posi- 
tion while this was still occulted, and data were recorded starting 
about thirty seconds before the nominal reapperance time. In re- 
ality, predictions resulted accurate to few seconds in all cases. 
Thanks to the work of Paranal software engineers, the event 
times were included in the files that control telescope and instru- 
ment operation (so-called Observing Blocks or OBs). These lat- 
ter were then loaded well in advance, and data acquisition would 
start automatically at the preset time. 

Two factors affected significantly the quality of the resulting 
data. Firstly, the active optics corrections of the primary mirror 
had to be switched off". As a result, the image quality was rather 
poor especially in the first part of the run. Towards the end of the 
run, we activated the active optics a few times in-between events, 
thus obtaining much better image quality. Secondly, in a few 
cases part of the signal fell outside the sub-window (clipping). 
Blind pointing was done by means of off'sets from a nearby star 
but, given also the bad image quality and the small detector sub- 
window, it was sometimes not sufficiently accurate. The sub- 
window was kept to 64x64 pixels, or w 9.5" squared, in order to 
achieve sufficiently fast read-outs. We stress that in general this 
does not have too many negative effects on the occultation light 
curve except of course for the signal-to-noise ratio (SNR). Fig.|3] 
illustrates the above mentioned eff'ects. 

The critical parameter for the angular resolution achieved in 
a LO, along with the signal-to-noise ratio (SNR), is the sampling 
of the light curve. For this, the burst mode implemented on the 
Aladdin detector of the ISAAC instrument was the critical factor, 
allowing us to obtain times of just 3.2 ms on a 32x32 subarray 
and 6.4 ms sampling on a 64x64 subarray. In order to increase 
the accuracy of source centering during blind pointing of the oc- 
culted sources, all the events but the first two were recorded in 
64x64 mode. The length of the data stream was kept between 
8000 and 14000 frames (see Table [TJ- We do not describe here 
the details of the burst mode for ISAAC, for which we refer to 
Richichi et al. (2008). The burst mode is also possible at several 
other ESO instruments both at Paranal and La Silla (see for ex- 
ample Domiciano de Souza et al. 120081 Poncelet et al. 120071 
Sicardy et al. I2007I I. The raw data need to be reformatted in 
proper FITS files, which are then processed using the AWLORP 
pipeline. Schematically, this latter performs both an automated 
light curve extraction based on a mask estimation which pre- 
serves the object-pixel connectivity, and computes an automated 
estimation of the basic parameters of each light curve, namely 
time of occultation, stellar signal, background signal, and rate of 
lunar motion. This is achieved via a novel algorithm based on 
wavelet-decomposition. Details on AWLORP can be found in 
Fors $20061 and Fors et al. (2008). 

Once the data are prepared with the proper structure, 
AWLORP runs automatically, saving the user from a lengthy and 
error-prone screening and producing preliminary fits and plots 
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Fig. 1. The area close to the sources occulted by the Moon in the night of March 21, 2006. The small dots represent sources in the 
2MASS catalogue with magnitude K< 7.5. The heavier dots are the sources for which we could record a lunar occultation. The 
apparent path of the Moon, moving from West to East, is also shown. The cross marks the Galactic Center. 
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Fig. 3. Illustration of the adverse effects of image clipping due to inaccurate blind pointing (left) and of poor image quality due 
to the active optics being switched off (middle, note that the target is a single star), as described in the text. The right figure is an 
example of a more regular case, for comparison. All three figures are time averages over several thousand frames on the subwindow 
of the Aladdin array in ISAAC. The windows are w 9.5" on the side. The 2MASS names of the occulted stars are indicated. 



for all Ught curves. A total of 30 events were found to have data 
of sufficient quality: the remaining 2 1 were discarded because of 
problems in the pointing or in the data acquisition. For the data 
analysis, we employed the model-dependent fitting ALOR (see 
Richichi et al. 1996) and the model-independent estimation of 
the brightness profile CAL (Richichi il989i 



3. Results 

The complete list of recorded occultation events is given in 
Table [T] It includes comments about the quality of the data and 
missed events. The K-band magnitudes range from 4. 1 to 7.4. In 
the table, and throughout the paper, we use 2MASS names. Most 
of the sources do not have a counterpart in the Simbad database. 
Only 14 of the sources reliably recorded can be cross-identified 
with optical or IR sources, and these are also listed in the table. 

We have tested the resolved/unresolved character of all the 
sources, using the method first outlined by Richichi et al. d 19961 1 
and applied in some of our previous papers such as in Fors et 
al. (2004) and Richichi et al. (2006a). The results are shown in 
Fig- 111 and we discuss below separately the resolved and the un- 
resolved sources. 

We show in the following just two examples of recorded light 
curves along with their best fits. The complete set of light curves 
and fits is available online. In spite of the previously mentioned 




SNR 

Fig. 4. Sources for which an upper limit on the angular diameter 
could be established using the method described in the text are 
marked with outlined circles. The filled circles are sources for 
which an angular diameter could be formally established, but 
for which the associated errorbar is such that they are undistin- 
guishable from being unresolved. The dashed line is an arbitrary 
representation of a tentative SNR-limiting angular resolution re- 
lationship. The triangles are sources for which the angular diam- 
eter and associated error lies tentatively or convincingly above 
the limiting angular resolution. 
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Table 1. List of the recorded occultation events 



2MASS id 


Frame 


NDIT 


Time 


K 


J-K 


Sp 


Cross-Id 


Notes 




(px) 




(UT) 


(mag) 








17383314-2836158 


32 


14000 


06:11:57 


5.64 


1.06 


- 


- 


df, 50% off edge 


17383141-2812445 


32 


14000 


06:15:19 


6.78 


2.09 


- 


- 


df 


17394834-2809136 


64 


12000 


06:23:51 


6.48 


2.33 


- 


- 


near terminator 


17385533-2836105 


64 


12000 


06:28:10 


7.34 


2.62 


- 


- 


df, 50% off edge 


17385407-2833227 


64 


12000 


06:32:56 


7.13 


2.67 


- 


- 


df, 50% off edge 


17391567-2826194 


64 


12000 


06:47:58 


6.46 


2.65 


- 


- 


df 


17394022-2820124 


64 


12000 


06:55:09 


7.12 


2.69 


- 


- 


df, close to edge 


17394602-2822089 


64 


12000 


07:00:05 


7.05 


2.67 


- 


- 


df 


17395038-2824061 


64 


12000 


07:04:00 


5.78 


1.10 


- 


CD-28 13409 


not seen 


17410318-2814151 


64 


12000 


07:08:12 


7.02 


2.57 


- 


- 


edge (75% off) 


17402920-2842361 


64 


12000 


07:14:29 


7.50 


2.73 


- 


- 


not seen 


17405112-2817130 


64 


12000 


07:18:35 


4.87 


2.20 


M6.5 


[RHI84] 10-333 


df 


17412116-2816010 


64 


12000 


07:22:54 


5.84 


1.55 


M5 


[RHI84] 10-356 


df, 2 stars 


17402750-2833099 


64 


12000 


07:26:35 


7.24 


2.67 


- 


- 


not seen 


17411415-2818051 


64 


12000 


07:29:49 


6.87 


2.00 


M5 


[RHI84] 10-352 


df, 2 stars 


17404347-2827185 


64 


12000 


07:33:08 


6.89 


1.98 


- 


- 


df, 50% off edge 


17410937-2822418 


64 


12000 


07:39:35 


6.92 


2.49 


- 


DENIS-P J174109.4-282242 


df 


17412349-2821538 


64 


12000 


07:44:18 


7.37 


2.39 


- 


- 


df 


17411109-2840486 


64 


12000 


07:47:31 


7.31 


3.06 


- 


- 


not seen 


17423294-2818473 


64 


12000 


07:52:21 


3.84 


1.27 


G8Iab 


HD 160706 


near terminator 


17414641-2822593 


64 


12000 


07:57:10 


5.09 


2.23 


- 


ISOGAL-P J174146.5-282301 


df 


17413435-2829225 


64 


12000 


08:01:38 


4.85 


2.11 


- 


DENIS-P J174134.4-282922 


df, 10% off edge 


17414185-2842287 


64 


12000 


08:05:22 


6.22 


2.66 


- 


- 


not seen 


17415470-2826596 


64 


12000 


08:09:09 


7.08 


2.54 


- 


DENIS-P J174154.7-282659 


df, 95% off edge 


17425620-2820370 


64 


12000 


08:13:32 


4.56 


1.39 


- 


2MASS J 17425620-2820370 


not seen 


17415719-2834236 


64 


12000 


08:17:55 


7.12 


2.89 


- 


- 


very faint 


17420509-2833465 


64 


12000 


08:22:10 


6.15 


2.88 


- 


- 


df, 50% off edge 


17423746-2825311 


64 


12000 


08:27:23 


5.20 


1.74 


- 


- 


df, very faint 


17424039-2826255 


64 


12000 


08:31:04 


4.56 


2.79 


- 


IRAS 17395-2825 


not seen 


17432345-2853503 


64 


12000 


08:38:21 


4.97 


3.46 


- 


ISOGAL-P J174323.6-285349 


not seen 


17432585-2823598 


64 


12000 


08:44:05 


5.92 


1.24 


- 


- 


df, 15% off edge 


17430791-2828016 


64 


12000 


08:48:53 


7.02 


3.55 


- 


- 


df 


17425497-2836593 


64 


12000 


08:53:18 


6.72 


3.77 


- 


- 


df, 50% off edge 


17432661-2827590 


64 


12000 


08:58:32 


6.26 


2.21 


M2: 


[RHI84] 10-443 


df, 50% off edge 


17430901-2837039 


64 


12000 


09:01:53 


6.99 


3.60 


- 


- 


df, 30% off edge 


17431348-2837004 


64 


12000 


09:04:34 


6.74 


3.72 


- 


- 


df 


17434830-2828129 


64 


12000 


09:10:32 


4.84 


2.29 


- 


- 


df, 95% off edge 


17435893-2827521 


64 


12000 


09:15:00 


5.97 


2.54 


- 


- 


df, second star 


17444369-2825027 


64 


12000 


09:17:24 


7.38 


3.21 


- 


- 


df, 5% off edge 


17435704-2830432 


64 


12000 


09:21:47 


5.61 


4.16 


- 


IRAS 17407-2829 


df, 25% off edge 


17434681-2840176 


64 


12000 


09:26:43 


4.45 


2.69 


M6.5 


ISOGAL-P J174346.6-284017 


not seen 


17435399-2841285 


64 


12000 


09:31:28 


4.11 


2.99 


- 


IRAS 17407-2840 


df 


17445251-2826479 


64 


12000 


09:36:45 


5.09 


3.10 


- 


- 


not seen 


17441366-2848501 


64 


8000 


09:39:53 


6.10 


0.99 


KO 


HD 316221 


5% off edge 


17443976-2830568 


64 


10000 


09:46:46 


7.40 


4.01 


- 


- 


not seen 


17452448-2826434 


64 


12000 


09:49:33 


6.72 


2.61 


- 


- 


not seen 


17450413-2853412 


64 


12000 


10:05:20 


5.75 


1.57 


M2 


[RHI84] 10-515 


df 


17452968-2829325 


64 


12000 


10:10:43 


6.64 


3.43 


- 


- 


gf, 2 more stars 


17445945-2843238 


64 


12000 


10:14:15 


4.97 


2.83 


- 


- 


gf 


17460150-2828035 


64 


12000 


10:19:24 


6.66 


4.47 


- 


- 


gf 


17453224-2833429 


64 


12000 


10:25:26 


5.31 


3.71 


- 


ISOGAL-P J174532.3-283338 


gf, 5% off edge 



Frame is the size of the subwindow and NDIT the number of frames used in the raw burst mode. The number of frames in the reformatted FITS cube is NDIT/2-I, and this also determines 
the total time span of the recorded data. In the notes, df stands for defocused and gf for good focus. 



problems, we note that our sample includes some of the best LO 
data ever recorded, both in terms of sensitivity and of SNR. The 
performance of the 8.2 m VLT telescope for LO work proved 
to be outstanding, not just due to its light collecting power but 
also to the great reduction of atmospheric scintillation - which is 
often the main limiting SNR factor for LO. In terms of angular 
resolution, the relationship shown in Fig.|4]is consistent with that 
established at other, smaller telescopes for which the detector 
integration time (DIT) was considerably shorter (see Richichi et 
al.[l996) ,2006aj . We note that the SNR of the fit shown in Fig.|5] 



is, as far as we know, the best ever achieved for a LO light curve. 
Just a few milliarcseconds away from the main star it would have 
been possible to detect a companion fainter than the primary by 
7.8 and 7. 1 mag in the occulted and unocculted part respectively, 
i.e. on opposite sides of the star This compares favorably with 
the best AO-assisted high-contrast imaging close to a bright star 
on this range of separations (Masciadri et al. |2005l l. 
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Table 2. Summary of results. 



(1) 
Source 



(2) 
|V| (m/ms) 



(3) 
VA',-1 



(4) 



(5) 
PA(°) 



(6) 

CAC) 



(7) 
SNR 



(8) 
Sep. (mas) 



(9) 
Br. Ratio 



(10) 
0UD (mas) 



17391567-2826194 


0.7702 


-2.9% 


-7 


268 


163 


211.0 


17412116-2816010 


0.3245 


5.7% 


-2 


350 


245 


105.4 


17435893-2827521 


0.4688 


4.9% 


-3 


321 


222 


233.2 


17444369-2825027 


0.1802 


62.3% 




see text 




68.2 


17435399-2841285 


0.6095 


-1.4% 


-5 


266 


168 


460.0 


17445945-2843238 


0.5945 


-1.7% 


-7 


263 


168 


251.4 


17453224-2833429 


0.5328 


7% 


-6 


302 


209 


126.9 



711.7 + 3.4 22.6 ±0.8 



171.6 ±25.8 47.8 ±7.5 



2.68 ± 0.66 

2.22 ± 0.37 

2.36 ±0.11 
1.91 ±0.30 

2.9 (shell) 



3.1. Resolved sources 

We have found two sub-arcsecond (projected) separation bina- 
ries, two stars with a convincingly resolved angular size, and 
three other stars which we consider as marginally resolved. 
These sources are listed in Table |2] using the same format al- 
ready used in Richichi et al. ( |2006al l and other papers in that 
series. In summary, the columns list the absolute value of the 
fitted linear rate of the event V, its deviation from the predicted 
rate Vt, the local lunar limb slope i//, the position and contact 
angles, the signal-to-noise ratio (SNR). For binary detections, 
the projected separation and the brightness ratio are given, while 
for resolved stars the angular diameter 0ud is reported, under the 
assumption of a uniform stellar disc. All angular quantities are 
computed from the fitted rate of the event. 

Several of our acquisition images show companions, but we 
report only two binaries which have projected separations below 
one arcsecond. This is rather large for the context of LO binary 
detections. In fact, for separations ^ (X'S the LO technique is rel- 
atively inaccurate, due to possible differences in the local limb 
slope at the points of occultation of the two components. The two 
objects are 2MASS 174121 16-2816010 and 2MASS 17444369- 
2825027. They are most likely projection and not physical pairs: 
we report them here for completeness, since no previous men- 
tion was present in the literature. The former was classified by 
Raharto et al. (1984) as a M5 star. We note that the difference be- 
tween predicted and observed lunar rate for the second star was 
quite significant, as reported in Table |2] The contact angle for 
this event was only 260°, i.e. nearly grazing, and the predicted 
lunar rate only -(X'l/s. Under these conditions, it is not uncom- 
mon that either a small error in the predictions or a local limb 
slope can significantly alter the computation of the effective PA 
and CA. 

The remaining stars in Table |2] are those which ap- 
pear significantly different from point-like. The two most in- 
teresting ones are 2MASS 17435399-2841285 and 2MASS 
17453224-2833429, which happen to be both masing AGB stars. 
Unfortunately the distance determination for these stars is uncer- 
tain, especially in the absence of detailed photometric and pul- 
sational properties. The former has an angular size which, albeit 
small in absolute terms, appears to be convincingly established 
(see Fig. |4] and Fig.|5]l. This star is cross-identified with IRAS 
17407-2840, a SiO maser previously measured by Messineo et 
al. ( 120021 120041 12005). Its radial velocity is negative and cannot 
be used for a kinematic distance, but Messineo (priv. comm.) 
estimates a distance of 3.7 kpc based on the average extinction 
and the model of Drimmel (12003). According to this, the angular 
size listed in Table |2] would translate to over 8 AU. Apart from 
consideration on the distance uncertainty, it is likely that the LO- 
derived size is not indicative of an uniform-disk diameter but 
rather a mix of the photosphere and the immediate circumstellar 
emission. 
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Fig. 5. Top: data (dots) and best fit (solid Hne) for 2MASS 
17435399-2841285. Bottom: enlarged by eight and displaced by 
arbitrary offsets for clarity, the residuals of the fits with a point- 
like source (above, x^ - 2.0) and with a uniform-disk model, 
resulting in a diameter of 2.36 + 0.11 mas (below, ;t'^ = 1.1). 



2MASS 17453224-2833429 has a similar natur e and set of 
references, with the addition of Sjouwerman et al. (2004). The 
fit for this source gave a SNR of 127 which is lower than ex- 
pected for the magnitude of this source (K=5 . 1 mag), and the UD 
size reported in Table |2] should be considered only as an indica- 
tion. A model-independent analysis by CAL revealed a faint but 
well defined extended emission around the central component, 
as shown in Fig. |6] The profile of this emission is almost sym- 
metrical, has the rim brightening typical of optically thin shells 
and is consistent with the masing nature of this star. In this case 
both a kinematic and an extinction distance can be estimated, at 
4.5 and 3.3 kpc respectively (Messineo, priv. comm.). By using 
their average, the inner shell radius of 10-15 mas would translate 
to ^40 AU. We note that the residuals of the uniform-disk model 
for 2MASS 17435399-2841285 indicate the possibiHty of a cir- 
cumstellar shell being detected also in this case (cfr. Fig.|5]l, but 
given that the x^ was already significantly good we did not pro- 
ceed with more extreme fits. The different signatures of the shell 
for these two maser stars could be related to their difference in 
J-K color (3.0 versus 3.7 mag). 

The remaining three stars in Table |2] appear resolved, al- 
though their location above the SNR-resolution relationship of 
Fig.|4]is less definite than for the two maser stars just mentioned 
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Fig. 6. Left: upper panel, data (dots) and best fit (solid line) for 
2MASS 17453224-2833429. The lower panel shows, on a scale 
enlarged by four and displaced by arbitrary offsets for clarity, the 
residuals of three different fits. Top, fit with a point-like source 
ix^ - 7.0). Middle, fit with a uniform-disk model, resulting in 
a diameter of 2.91 + 0.27 mas {x^ - 6.3). Bottom, fit with the 
model-independent CAL method (x^ - 1.6), by far the most 
convincing. Right: brightness profile reconstructed by the CAL 
method, corresponding to the residuals of the lowermost fit on 
the left. 
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Fig. 7. Plot of the SNR achieved in our fits of both resolved and 
unresolved source, as a function of the total measured flux. 



and we prefer to consider their diameters as tentative for the mo- 
ment. There are no known cross-identifications or Uterature en- 
tries for these stars. 



3.2. Unresolved stars and performance 

Most of the stars in our sample appear unresolved, as is to be 
expected from random observations of field stars. Richichi et 
al. (2006a) have discussed the incidence of binary stars in such 
samples for LO observations based on the 2MASS catalogue, 
and our present results are consistent within the uncertainty of 
small statistics. The unresolved stars are also useful, namely in 
that the relationship of Fig.|4]establishes rather strict upper limits 
on their size and on the brightness ratio of unseen components, 
if any. Thus these stars can be safely used as calibrators for long- 
baseline interferometry. 

Further, the unresolved stars allow us to investigate the per- 
formance of the burst mode at the VLT for LO, similar to what 
was done for other instrument and telescope configurations in 
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Fig. 8. Plot of the total measured flux against the K magnitude of 
the sources (from 2MASS). The line shows the expected counts, 
based on the ISAAC Exposure Time Calculator. Note the dis- 
crepant cases, due to the problems mentioned in the text. 



the past (see Fors et al.'2004). In Fig.|2]we plot the SNR achieved 
in the fits of all sources (the unresolved ones being the majority), 
against the total measured flux. This is a better quantity to use 
than the actual magnitude, since as mentioned earlier many of 
our sources were recorded at the edge of the subarray field with 
considerable loss of signal. Also intrinsic variability might be a 
cause for some of the stars. Indeed, Fig. |8] shows many outliers 
from the ideal K magnitude-observed flux are present. Using the 
measures signal instead, Fig.|2]shows a rather linear relationship 
across over two orders of magnitudes in flux. This implies that 
at the faint end we were still in a photon-limited regime, and at 
at the bright end the incidence of scintillation and possible non- 
linearities was not yet significant. We estimate that non-linearity 
of the detector should become important above Ka!3 mag, i.e. 
one magnitude above our brightest target for this run. 

In spite of the outlyers, the relationship between K magni- 
tude and counts clearly emerges from Fig. [8] and it is in excellent 
agreement with the expected ISAAC performance. Comparing 
the two figures, we estimate at first approximation that under 
the same conditions of sampling and DIT (6.4 ms) it should 
have been possible to obtain SNR=1 for Ka!l2mag. This is the 
minimum SNR necessary to detect a faint companion in a LO 
light curve, although measurements of angular diameters or ex- 
tended emission would require SNRSIO. We note that, at the 
faint magnitude end, critical parameters are the diffuse scattered 
light background, which is a function of the lunar phase and dis- 
tance of the occultaton point from the terminator, and the pixel 
scale. Ideally, this latter should match the seeing disc, but in our 
case with a pixel scale of about 0'.'148 we significantly oversam- 
pled the background. 

4. Conclusions 

We have presented the first LO results obtained using the so- 
called burst mode of the ISAAC instrument at the ESO Antu 
VLT telescope. This mode permits to record data streams on 
a small subarray at high temporal frequency (3.2 ms, although 
mostly 6.4 ms were used for the present work). We recorded 
5 1 occultation events over about four hours during a close ap- 
proach of the Moon to the Galactic Center The events were 
reappearances from behind the dark limb. This, coupled with 
initial problems of focussing and image quality due mostly to 
the commissioning nature of the run, resulted in several sources 
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being missed or without sufficient quality. Nevertheless, 30 of 
the events resulted in light curves of good, or even excellent 
quality, including some of the best SNR traces ever recorded 
with this technique. We have revealed two binaries, three stars 
with a marginally resolved angular diameter of the order of 
2 mas, and two resolved masing AGB stars which appear to be 
in the foreground of the galactic bulge. For one of these, 2MASS 
17453224-2833429, we were able to recover the brightness pro- 
file of the extended emission. This has the typical signature of 
an optically thin shell, with inner radius of order 10-15 mas, cor- 
responding to about 40 AU. Follow-up studies are required to 
obtain further information on these stars, given that they are all 
heavily reddened and most have no known counterpart or litera- 
ture entries. 

Our analysis, including the unresolved sources, has estab- 
lished an unprecedented performance of the LO technique in the 
near-IR using a very large telescope and a fast read-out mode. 
We find that the angular resolution varies between 4 and 1 mas, 
and even less in the case of very high SNR. Data quality is gen- 
erally impressive, compared to smaller telescopes, mostly due to 
the scintillation reduction effect of the large mirror. The dynamic 
range is also significantly improved with respect to previous ob- 
servations, and we have shown that in some cases it would have 
been possible to detect a companion fainter than the primary by 
about 8 magnitudes in broad K band on angular scales compara- 
ble to the Airy disk of the telescope. This result compares favor- 
ably with other AO-assisted high contrast imaging techniques. 

The random selection of sources to be occulted and the fixed- 
time nature of the events remain the main downsides of the LO 
technique. However, thanks to the very short telescope time re- 
quired for each observation, LO are an ideal filler program for 
times when the telescope is idle or other programs cannot be 
conveniently executed. The service mode observations imple- 
mented at the VLT provide the right context for this, and we have 
obtained approval for such a filler program during ESO observ- 
ing periods 80 and 8 1 (from October 2007 through September 
2008). Several hundreds of OBs have been prepared for each 
period, and are in the queue awaiting execution whenever five 
minutes of telescope time become available. 
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